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Thermal emission from WASP-33b, the hottest known 
planet* 
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ABSTRACT 

We report ground-based observations at 0.91 fim of the occupation of the hot 
Jupiter WASP-33b by its A5 host star. We measure the planet to be 0.109 ±0.030 per 
cent as bright as its host star at 0.91 /im. This corresponds to a brightness temperature, 
Tb — 3620^250 K, significantly higher than the zero-albedo equilibrium temperature 
for both isotropic re-radiation (2750 ± 37 K) and uniform day-side only re-radiation 
(3271 ±44 K), but consistent with the zero-redistribution temperature (3515 ± 47 K). 
This indicates that the heat redistribution from the day-side of WASP-33b to the 
night side is inefficient, and further suggests that there is immediate re-radiation, and 
therefore little or no redistribution, of heat within the day-side. We also detected the 
stellar pulsations of WASP-33, which we model as the sum of four sinusoids, with 
periods of between 42 and 77 minutes and amplitudes of 0.5 to 1.5 mmag. 
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1 INTRODUCTION 

Extra-solar planets which transit their host stars are also 
likely to exhibit occultations (secondary eclipses), in which 
the planet passes out of sight behind the star. Observing 
an occultation allows measurement of the planet's emergent 
flux. This flux consists of thermal emission and reflected 
light, but is strongly dominated by thermal emission in the 
infrared wavelength regime. The first detections of thermal 
emission from an exoplanet were made us i ng the Spitzer 
space telescope (|Charbonneau et al ] 120051 : iDeming et alj 
2005), since when the occultations of many more planets 
have been observed at wavelengths of 3.6 /im and greater 
with Spitzer. 

Complementary to Spitzer observations at 3.6, 4.5, 5.8, 
8 and 24 fim, are ground-based observations made in the 
near infrared. Such observations are important because they 
extend the observed planetary spectral energy distribution 
(SED) to lower wavelengths, and to the peak of the SED, 
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which is around 1 /im for the shortest-period, hottest plan- 
ets. 

Measurement of a planet's occultation depth allows the 
brightness temperature to be determined, which leads to an 
evaluation of how efficiently heat is redistributed from the 
day-side to the night-side of the planet. 

A planet's observed SED (from occultation observa- 
tions at several different wavelengths) allows the chem- 
ical compos i tion of the atmosphere to be inferred (e.g. 
I Swain et al.l |2010T ). Detection of a thermal inversion in 
the atmosphere is also possible; observations to date indi- 
cate two classes of hot- Jupiter atmosphere, those with and 
those without such an inversion. One possible explanation 
for the existence of these tw o class es is that provided by 
iKnutson. Howard &: Isaacson! l|2010l ) who argue that chro- 
mospherically active stars host planets with no temperature 
inversion. 

Occultation observations also allow refinement of the 
orbital parameters of a planetary system, particularly the 
orbital eccentricity, which is related to both the timing of 
the occultation and its duration. 

WASP-33 (HD 15082 ) was first iden t ified a s a tran- 
siting planet candidate by IChristian et al.1 (|2006l ), but was 
identified as fast rotator, which rules out the use of pre- 
cision radial velocity measurements to confirm its plane- 
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tary nature and to he l p char acterise the system. Instead, 
ICollier Cameron et al.l £2010), hereafter CIO, used line- 
profile tomography during transit to confirm WASP-33b as 
the first known planet to orbit an A-type star. 

Given the nature of its host star and its very short or- 
bital period (1.22 d), WASP-33b has the largest equilib- 
rium temperature (Teql ~ 2700 K, assuming zero-albedo and 
uniform heat redistribution; CIO) of any known exoplanet. 
WASP-33b is an excellent target for ground-based occulta- 
tion observations, due to the large predicted planet-to-star 
flux ratio. Planets with similarly large b rightness tempera - 
tures include WASP-19b (T cq i « 2000 K; lHebb et aDl2010h . 
occultations of which h ave been observed from the ground in 
-ff-ba nd and if-band (jAnderson et al.ll2010al; Gibson et al.l 
120101 ) and WASP-12b (T eq i « 2500 K; iHebb et al.ll2009h . 
observed from th e ground in the Ks~, H- and J-bands by 
ICroll et al.1 l|201ll ). 

C10 also reported that the stellar line profiles show ev- 
idence for non-radial pulsations, and suggest that the star 
may be a 7 Dor-type variable. Further, ph otometric, ev- 
idence for these pulsations is presented in iHerrero et al.l 
i|201ll) . hereafter Hll, who suggest that WASP-33 is a 
5 Scuti-type variable. 

In this paper, we present ground-based observations of 
the occultation of WASP-33b at 0.9/rni. 



2 OBSERVATIONS AND DATA REDUCTION 

We observed the occultation of WASP-33b on 2010 October 
29/30 using the Auxiliary Camera (AC AM) of the William 
Herschel Telescope (WHT). The AC AM CCD covers a field- 
of-view of 8' x 8' with 2048 x 2048 pixels, giving a plate scale 
of 0.25" per pixel. The target was positioned on the CCD 
such that three relatively bright comparison stars are present 
in the images. A total of 203 images were taken from 22:43 
to 05:55 UT. 

The observations were conducted using a narrow-band 
S[III] filter (ING filter #S9077), which has a central wave- 
length of 9077 A, and a FWHM of 54 A. The telescope was 
defocussed to spread the light from the target over a large 
number of pixels, thus reducing noise associated with im- 
perfect flat-fielding and allowing longer integration times. 
Exposure times of 100, 120 and 150 s were used. The con- 
ditions were not photometric; thin cirrus cloud affected the 
observations for the duration of the night. 

Aperture photometry was performed on the flat-field 
a nd bias corrected images , using the photometry pipeline 
of ISouthworth et al.l (|2009l ). which is based around the AS- 
trolib / APER routinqj. The pointing was monitored by 
cross-correlating each image with a reference image, and the 
apertures shifted to follow the stars. 

Photometry was performed using a range of different 
aperture radii, from 8 to 65 pixels. A radius of 50 pixels was 
chosen to maximise the target signal-to noise, and ensure 
that a large fraction (> 95 per cent) of the target flux is 
included in the aperture. The flux from the three compari- 
son stars was combined to form an ensemble reference star, 
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Figure 1. Normalised raw light curve. 



approximately 0.5 magnitudes fainter than WASP-33. The 
resulting differential light curve is shown in Fig. [T] 



3 DATA ANALYSIS 

In addition to an apparent eclipse-like signal (centred about 
JD ~ 2455499.58 as expected, assuming a circular orbit), 
the raw light curve (Fig. [TJ exhibits large systematic ef- 
fects, specifically (i) a series of sine-like oscillations, which 
we attribute to the stellar pulsations for which evidence was 
presented by C10 and Hll; and (ii) an increase in flux across 
the second half of the light curve. 

We analyse the new occultation data alongside the 
SuperWASP, Keele and JGT light curves and the ra- 
dial velocities of C10, and the Montsec Observatory light 
curve of Hll. This global analysis is performed using 
th e Markov Chain Monte Carlo (MCMC) code desc ribed 
in ICollier Cameron etafl \200H ), Enoch et all (|2010h and 
lAnderson et al.1 |2011bT v 

In addition to fitting the occultation, we also simul- 
taneously fit functions describing the stellar pulsations and 
systematic effects present in the raw light curve (Fig.[T]). The 
de-tren ding within the MC MC code is done using the SVDFIT 
routine (jPress et al.ll 19921 ). Initially, we neglected the stellar 
pulsations, and attempted only to remove the instrumen- 
tal effects. We tried de-trending with several parameters: 
airmass, time, sky background amplitude, and the position 
of the target on the CCD; we also tried de-trending using 
various combinations of parameters and various functional 
forms for the fit to each parameter. We used the Bayesian 
information criterion (BIC) to discriminate between the dif- 
ferent models. In determining the best de-trending model we 
did not re-scale the error bars on any of the data. The low- 
est BIC value was produced when the data were de-trended 
with a quadratic function of the sky background, 



g — ao + ai/ s ky + a2/ s ky 



(1) 



further details see http://idlastro.gsfc.nasa.gov 



where / s k y is the sky background flux and ao, a\, and 02 are 
the coefficients we fit for. 

The light curve resulting from our de-trending with sky 
background, while an improvement upon the raw light curve, 
still exhibits the sine-like variations present in the raw light 
curve which we attribute to stellar pulsations. It is apparent 
the signature of these pulsations is multi-periodic; a model 
consisting of a single sine curve produced a very poor fit to 
the data. 
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Figure 2. Periodogram output from Period04, for the out-of- 
occultation data. 

We conducted a Fourier analysis on the out-of- 
occultation data using the Period04 time-series analysis 
software package (|Lenz fc Bregerl feOQS). This reveals three 
frequencies which are deemed significant (with a signal-to- 
noise ratio greater than 4). These frequencies correspond to 
periods of 53.99±0.29 min, 76.90±0.62 min and 41.85±0.31 
min (see periodogram, Fig. [2}. The 68-min period of Hll is 
present in our periodogram, but with a slightly lower ampli- 
tude than the 54-min period. When prewhitening with the 
54-min period, the the significance of the 68-min period is 
much reduced, suggesting that one of these two periods may 
be an alias of the other. When, instead we used the Hll 
period and amplitude as the starting point for a Period04 
analysis yields a different set of frequencies, but a poorer 
overall fit to our data. We also tried a period search of our 
data combined with the existing photometry, but the results 
from this were inconclusive. 

We then fitted for either a single pulsation period (that 
of Hll), three pulsation periods (from our Period04 anal- 
ysis) or four pulsations periods (our periods and that of 
Hll) with an MCMC code (see below). Comparing the 
BIC values of the resulting models we found that the 4- 
period model (ABIC = 0) is preferred to both the 3-period 
(ABIC = 391) and the 1-period (ABIC = 1057) models. 

The four periods and uncertainties were input to a ver- 
sion of our MCMC code, modified to include the four pul- 
sation periods as proposal parameters. The periods were, 
therefore, allowed to vary but Gaussian priors were applied 
to prevent the periods from drifting too far from the fitted 
periods. The priors were applied by means of a Bayesian 
penalty added to our merit function (x 2 ), given by 

BP Pi = (Pl ~ Plfi)2 (2) 
a P,,o 

where P^o and <xf>. are the initial value and uncertainty 
of a particular period, Pi, respectively. The amplitudes and 
phases of each of these sine waves were fitted using SVDFIT, 
along with the coefficients of Equation[T] The function fitted 
is of the form 




where g(f s ky) is Equation [T] t is time, and bi and d are the 
coefficients we determine in the fit. This is equivalent to 

h = Si/sky) +J2 A ' Shl + Si ) W 



Table 1. Comparison of de-trending models, g, (see Equation [3). 
ABIC is the value of BIC relative to that of the model with the 
lowest value of BIC 



Functional form of g ABIC 



ao - 


l-ai/sky + a2/ s 2 ky 





ao - 


h <H/sky 


1 


a - 


h a 1 St + a 2 (6t) 2 


11 


ao - 


^l/sky 


12.3 


ao - 


- ai sec z + a2 sec 2 z 


118 



/sky is the sky background flux, z is the zenith distance of the 
target, and St is the time since the start of the observations. 

where Ai are the amplitudes, and 5i the phase offsets of each 
sine term. 

We then again tried different forms of Equation [T]to de- 
trend the data, again finding that a quadratic in sky back- 
ground results in the lowest BIC value (Table [TJ. Once we 
had selected the functional form of our trend model, we per- 
formed a further MCMC analysis, this time scaling the error 
bars of each photometric dataset so as to obtain a reduced 
X 2 of unity. 

We initially fitted for the orbital eccentricity and the re- 
sulting solution is significantly eccentric (e = 0.174to'o3g), 
but ecosu is close to, and consistent with, zero (ecosw = 
0.0027^ 

^0033)- The argument of periastron, w, in this solu- 
tion is — 89.0if'2°, i.e. the major axis of the orbit is aligned 
almost exactly with our line-of-sight. We suggest that this is 
an improbable configuration, and is caused by a degeneracy 
in fitting the occultation and the stellar pulsations, ecosw 
is constrained by the timing of the occultation, whereas the 
duration of the occultation informs us about esinw. 

Furthermore, in the case of WASP-33, the shape of the 
radial velocity curve places almost no constraint on the ec- 
centricity (because the RVs are imprecise due to the nature 
of the star). This allows the fitted occultation duration to 
deviate from the transit duration to a large extent. The light 
curve, after de-trending for sky background and stellar pul- 
sations still clearly contains correlated ('red') noise, which 
is likely to be caused by pulsation activity at further peri- 
ods (Fig. [4]). It is this red noise that mimics the ingress and 
egress of the occultation and results in an occultation dura- 
tion significantly shorter than the transit duration, and an 
argument of periastron close to —90°. 

The eccentric solution requires that the planet be mov- 
ing more slowly during transit than for the circular solu- 
tion, and this in turn requires that the stellar radius is 
smaller in the eccentric case. This results in a stellar density 
(p* = 0.75^g'}2 P©) which is greater than that of a main- 
sequence star with the mass of WASP-33 (p* w 0.5 P0; Gray 
2008). The stellar density of our circular solution, however, 
lies very close to that prediction - p* = 0.44 ± 0.05 Pq. 

Although short-period planets have been found in sig- 
nificantly eccentric orbits (e.g. WASP-14b. |joshi et al . 2009; 
iHusnoo'etaflteOld ). the orbit of WASP-33b is expected to 
have been circularised by tidal interactions with the star. 
These are anticipated to be particularly strong because of 
the relatively large radii of planet and star, and the small 
orbital separation. 

Using Equation (1) of Ijackson. Greenberg fc Barnes! 
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Table 2. System parameters for WASP-33 



Parameter (Unit) 



Orbital period, P (d) 

Epoch of mid-transit, T c (HJD) 

Transit and occultation durations, T14 = Tsg (d) 

Transit and occultation ingress and egress 

durations, T 12 = T34 = T 56 = T 78 (d) 
Transit depth, AF = Rp /R2 
Impact parameter, b 
Orbital inclination, i (°) 

Orbital eccentricity, e 

Occultation depth, AFo.gi^m 

Stellar effective temperature, T„ e ff (K) 

Stellar metallicity, [M/H] 

Stellar mass, M, (M ) 

Stellar radius, ii* (-Rq) 

Stellar surface gravity, log g„ (cgs) 

Stellar density, p* (po) 

Planet mass, Mp (Afj up ) 

Planet radius, Rp (i?j up ) 

Orbital major semi-axis, a (AU) 

Planet equilibrium temperature, T eq i(A_o f = i) (K) 

Planet equilibrium temperature, T eq i(^_g f_2) (K) 



Planet equilibrium temperature, T 



eql(A = 0,f=|) 

Planet brightness temperature, Tp^o.gipm (K) 
Pulsation period 1, Ppi (min) 
Pulsation period 2, Pp2 (min) 
Pulsation period 3, Pp3 (min) 
Pulsation period 4, Pp4 (min) 
Pulsation amplitude 1, Api (mmag) 
Pulsation amplitude 2, Ap2 (mmag) 
Pulsation amplitude 3, Ap3 (mmag) 
Pulsation amplitude 4, Ap4 (mmag) 



(K) 



(CIO) 



This work 



1.2198669 ±0.0000012 
2454163.22373 ± 0.00026 



1.21986983 ± 0.00000045 
2454590.17948 ± 0.00028 
0.11224 ±0.00084 



0.01136 ±0.00019 
16+ 010 
87.7 ± 1.8 



7430 ± 100 

0.1 ±0.2 
1.495 ±0.031 
1.444 ±0.034 

4.3 ±0.2 
0.497 ±0.024 

< 4.1 
1.497 ±0.045 
0.02555 ±0.00017 



01149+ (u)00!,7 

U.Ullia_ .()|)()S4 

01041+ 00023 
U.U1U11_ .00021 

0.32j 

84.9^ 



,+0.09 
-0.12 
,+2.1 
-1.7 



(fixed) 
0.00109 ± 0.00030 
7430 ± 100 
0.1 ± 0.2 
1.512 ± 0.040 



1 cri 9 +0.060 
1.0iz_o 054 
a+0. 027 
-0.029 



4.258Z 
0.437 ± 0.045 



< 4.59 (3-cr) 
1.501 



+0.073 



-0.064 

0.02565 ± 0.00023 
2750 ± 37 
3271 ± 44 
3515 ± 47 

F200 



36201 250 
« 5+09 

,+2.6 
-4.0 



76.3Z 
41.9 ±0.5 

1.479 ±0^069 
0.567 ±0.134 
0.766 ±0.115 
0.605 ±0.105 



(2008), we calculate the circularisation time-scale, r e = 
(ell) using the values of M„, R,, a, and Rp quoted in 
Table [2] and Alp = 4.1Mj up (the upper limit established 
by C10, which will maximise the calculated circularisation 
time-scale). We obtain 

( 0.557 0.024 \ 

where Qp and Q* are the tidal dissipation parameters for 
the planet and star respectively. Adopting Qp = 10 5 ' 5 
and Q t = 10 6 5 (the best-fitting v alues from the study of 
I Jackson. Greenberg fc Barnedliool ). we find r e = 1.72 Myr. 
This time-scale is an order of magnitude less than the 25 
Myr main-sequence age estimate of C10, and a very small 
fraction of the 500 Myr upper limit to the age. Although 
this suggests that there has been ample time for the orbit to 
have been circularised, r e is affected by large uncertainties 
in the Q parameters, so we cannot say for certain that this 
is the case. 

C10 determined that the orbital axis of WASP-33b is 
misaligned with respect to the stellar spin axis, such that the 
planet orbits in a retrograde sense. This suggests that the 
system may not exist in a dynamically relaxed state, and the 
orbit may be eccentric. The time-scale for co-planarisation 



of the orbit, however, is very much long er than that for its 
circularisation. Using Equation (14) of I Winn et alj (|2005l ) 
we obtain the characteristic time-scale, r^,, for significant 
change in the angle ip between the planetary orbital angular 
momentum vector and the stellar rotation angular momen- 
tum vector, 

^ = 1L4 (l^) Gyr ' (6) 

Here, we have assumed values of 0.01 for the apsidal motion 
constant, V0.1 for the dimensionless radius of gyration, and 
a planet mass of 4.1 Mj up . The time-scale increases with 
decreasing planet mass. It is therefore perfectly reasonable 
to suppose that the orbit of WASP-33b is circular, despite 
the fact that it is retrograde. 

Furthermore, the empirical evidence provided by other 
planetary system s lends support to the circular hypothesis. 
I Pont et al.l (|201ll ) present an analysis of the eccentricities of 
a large number of hot Jupiters, determining that the only 
known planets with significantly detected, large eccentrici- 
ties are massive and have relatively long-period orbits. Sev- 
eral planets with large spin-orbit misalignments are known 
to have near-circul ar orbits, such as WAS P-17b which has 
a retrograde orbit (| Anderson et al.ll2010bl ). and an orbital 
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Figure 3. Stellar pulsations. Photometry, with sky background 
and occultation models subtracted (blue points), overplotted with 
our best-fitting pulsation model (red curve). 
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Figure 4. Occultation photometry, de-trended for sky back- 
ground and with stellar pulsations subtracted, normalised to 
unity in occultation (blue points), with best-fitting model over- 
plotted (red curve). 

eccent ricity close to zero, e = 0.028to'ois <|Anderson et al.l 
l2011ch . 

All of the above evidence leads us to rule out the eccen- 
tric solution, and to t ake the approach (advocated by e.g. 
lAnderson et al.ll2011al ) of adopting a circular solution. The 
best-fitting parameters for this circular model are presented 
in Table [2] The fit to the pulsations is shown in Fig. [3] and 
the fit to the occultation is shown in Fig. [3] 

We assessed the presence of correlated noise in the oc- 
cultation light curve residuals, by plotting the RMS of the 
binned residuals to the best-fitting model as a function of 
bin width (Fig. [5]). The RMS of the binned residuals devi- 
ates significantly from the uncorrelated noise expectation, 
indicating significant levels of red noise are present in the 
light curve. This is likely to be due to the presence of pulsa- 
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Figure 5. RMS of the binned residuals as a function of bin size 
(solid, blue curve). The dashed, red curve is the expected RMS 
for noise that is purely uncorrelated, corresponding to the RMS 
of the unbinned data divided by the square root of the number 
of points per bin. 



tion signals in addition to those we fitted. The results of a 
Fourier analysis to determine whether any further periodic- 
ity is present in the data is inconclusive, in part due to the 
limited baseline of the data. To determine whether fitting an 
occultation model to the data is justified, we performed an 
F-test to compare our best-fitting occultation model with a 
constant-flux model. The resulting p = 0.013 (i.e. there is 
only a ~ 1% chance of our occultation fit arising if there 
were no underlying occultation) indicates that the occulta- 
tion model is preferred at the > 95 per cent level. 

The remaining correlated noise present in the occul- 
tation light curve means that the uncertainty on the oc- 
cultation depth is likely to have been underestimated. To 
calculate a realistic uncertainty, we emp loy the residual per - 
mutation or 'prayer bead' method (e.g. iGillon et all l2007l l . 
We took the light curve residuals and shifted each resid- 
ual to the subsequent observation before adding back the 
trend, pulsation and occultation models. This process was 
repeated, shifting each residual two observations, and so on, 
resulting in a total of n b s — 1 = 202 light curves, each of 
which have the time-structure of the correlated noise pre- 
served. A separate MCMC analysis was performed on each 
of these light curves. 

We then calculated the median of the best-fitting oc- 
cultation depths and pulsation periods from the new fits, as 
well as the uncertainty limits which enclose 68 per cent of 
the values around the median. The uncertainties on these 
parameters quoted in Table [2] are those calculated from the 
permutation of the residuals. In all cases, the median values 
from the residual permutation analysis are consistent with, 
and close to the best-fitting values from our single MCMC 
analysis. As expected, the uncertainties on these parameters 
derived from the permutation analysis are larger. The origi- 
nal occultation depth we found was 0.00108±0.016, whereas 
the value from the residual permutation is 0.00109 ± 0.030. 



4 DISCUSSION 

To calculate the brightness temperature of the planet, we 
modelled the star using a synthetic spectrum of an A5 star 
I Pickles! (1 19981 ) normalised to the integrated flux of a 7430 K 
(Table [2| black-body emitter, and the planet as a black- 
body of temperature Tb. We then defined the measured oc- 
cultation depth as the product of the ratio of the bandpass- 
integrated plane^Bd_s_ta£ihrxes_md the planet-to-star area 
ratio (e.g. ICharbonneau et al.ll2005l ). We calculate a bright- 
ness temperature, Tb = 36201250 K. The uncertainty on Tb 
only accounts for the uncertainty on the measured occul- 
tation depth, which is the dominant source of error. Addi- 
tional, smaller sources of uncertainty on the ratio of plane- 
tary and stellar areas, and the stellar effective temperature 
are not accounted for. 

The equilibrium temperature for a zero-albedo planet 

is given by T p ,a=o = T*,eff \f§^, where / = 1 in- 
dicates isotropic re-radiation over the whole planet (i.e. 
the redistribution of heat from the day-side to the night- 
side is fully efficient). The case where the heat is uni- 
formly distributed on the day-side of the planet, but there 
is no redistribution of heat to the night-side corresponds 
to / = 2. A third case, corresponding to / = f (equiva- 
lent to e = and / = | in the notations of I Cowan &: Ago] 
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l201ll and lL6pe z- Morales fe Seager]|2007l respectively), occurs 
when the incident radiation is immediately re-radiated and 
so there is no redistribution of heat even within the day-side. 
Because the hottest region of the day-side is most visible 
close to occultation, a deeper occultation is observed than 
if the day-side had a uniform temperature. 

The brightness temperature of 3620 K is greater than 
either the equilibrium temperature for uniform or no redis- 
tribution to the night-side (respectively 2750 K and 3271 K, 
Table |2J. The brightness temperature is consistent with the 
zero-redistribution (/ = |) case (3515 K), however. This 
suggests that the heat-transport from the day-side to the 
night-side is inefficient, and further that the day-side is un- 
likely to have a uniform temperature. This is in line with 
observations of other highly irradiated planets, which are 
also found to have p oor heat redistribution efficiencies (e.g. 
ICowan fe Agofolll ). 



5 CONCLUSIONS 

We have detected thermal emission from WASP-33b. We 
measure the occultation depth at 0.91/im to be 0.109±0.030 
per cent, which corresponds to a brightness temperature of 
3620^250 K, the hottest such temperature recorded for an 
exoplanet. 

We also detect the non-radial pulsations of the host star 
in the photometry; and conclude that a multi-periodic solu- 
tion is required to fit the pulsation signal. Our best-fitting 
model of the pulsation signal consists of sine terms with pe- 
riods of 53.5, 76.3, 41.9, and 66.6 min, with amplitudes of 
1.5, 0.6, 0.8 and 0.6 mmag respectively. None of these pe- 
riods corresponds directly to the period of 68 min fitted by 
Hll, but are similar, and our greatest amplitude agrees with 
their 0.9 mmag. We do however find evidence for the Hll 
68-min period in our periodogram, and it appears that the 
54-min and 68-min periods may be aliases of each other. We 
do not claim that these periods are definitive, but rather 
they are the best fit to our data. Our analysis suggests that 
the pulsations of WASP-33 are complex and multi-periodic 
in nature. This multi-periodic, ~ 1 h, pulsation signature 
lends support to the conclusion of Hll that WASP-33 is a 
5 Scuti-type variable. Our ability to determine the true pe- 
riodicity of the pulsations is limited by the relatively short 
baseline of our data. 

Although fitting for the orbital eccentricity returns a 
significantly non-zero value of e, we argue that because 
ecosw is essentially zero, the eccentric solution is improba- 
ble. Furthermore, the eccentric solution is incompatible with 
the stellar analysis of C10, so we therefore adopt a circular 
model for the orbit. 

High signal-to-noise photometry with a longer base- 
line is required to study better the pulsations, to determine 
the complete nature of the periodicity and to investigate 
whether, as Hll suggest, there is evidence for star-planet 
interactions in the WASP-33 system. Such a characterisa- 
tion of the pulsation periods would also allow the pulsation 
signal to be more cleanly subtracted from occultation pho- 
tometry. This, in turn, would allow a more precise measure- 
ment of the occultation depth and the unresolved questions 
concerning the orbital eccentricity to be answered. Measure- 
ments of the occultation depth at different wavelengths to 



the one presented here are required in order to construct a 
spectral energy distribution and hence to characterise the 
planetary atmosphere. 
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